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^ ! Abstract 



m 
> 

psj ' parameter field is derived for the case of point defects in the 0{n) symmetric 



An explicit expression for the vortex velocity field as a function of the order 



time-dependent Ginzburg-Landau model. This expression is used to find the 

' vortex velocity probability distribution in the gaussian closure approximation 

in the case of phase ordering kinetics for a nonconserved order parameter. 

The velocity scales as L~ l in scaling regime where L ~ i 1 / 2 and t is the time 

q ' after the quench. 

O ■ 

The importance of the role of defects in understanding a variety of problems in physics 
is clear. In certain cosmological |L] and phase ordering problems key questions involve 



an understanding of the evolution and correlation among defects like vortices, monopoles, 
disclinations, etc. In studying such objects in a field theory questions arises as to how one 
can define quantities like the density of vortices and an associated vortex velocity field. 
The purpose of this paper is to identify the appropriate vortex- velocity field in the context 
of an 0(ri) symmetric time-dependent Ginzburg-Landau (TDGL) model for the case of 
point defects where n = d and d is the spatial dimensionality. Using this rather general 
definition for the velocity field the distribution of velocities is determined in the case of the 
late state phase ordering using the gaussian closure approximation for a nonconserved order 
parameter. The physical results are that the velocity scales as Lit)' 1 where L(t) ~ t 1 ^ 2 is 

1 



the characteristic scaling length for the order parameter correlation function which grows 
with time t after the quench. The vortex velocity probability distribution function is given 
in this approximation by 

p(| + 1) l 

P ( W o) = ( n y2\n/2 7 \ (n+2)/2 W 

( 1 + {v Q f/vA 

where the parameter v is defined below and varies as L~ l for long times. 

The focus here is on the defect dynamics generated by the TDGL model satisfied by a 
nonconserved n-component vector order parameter tp{f,t): 

where T is a kinetic coefficient, F is a Ginzburg-Landau effective free energy assumed to be 
of the form 

d d r{^{V^) 2 + Vm) (3) 

where c > and the potential is assumed to be of the degenerate double-well form, if is a 
thermal noise which is related to T by a fluctuation-dissipation theorem. 

Consider a system with n = d where there are topologically stable point defects || 
formed in a phase ordering system (quenched for example from a high temperature disordered 
state to a temperature below the order temperature). As pointed out by Halperin [|J, and 
exploited by Liu and Mazenko ||, the vortex density for such a system can be written as 

P =mv (4) 

where T> is the Jacobian (determinant) for the change of variables from the set of vortex 
positions rj(t) (where ip vanishes) to the field ip: 

where e W)A t 2 ,...,At„ is the n-dimensional fully anti-symmetric tensor and summation over re- 
peated indices here and below is implied. 



The first goal here is to derive the equation of motion satisfied by p. Toward this end 
one needs two identities whose proof is relatively straightforward. The first identity is given 
by: 

dT> 

— = V a Ji K) Identity I (6) 
where, for a general vector A, the current jf^ is defined as 

Identity I is just a statement that the determinant D is a conserved invariant. Notice that 
the superscript K on J in this identity is defined by the right hand side of Eq.(H). The 
second identity takes the form for general vector A: 

J { a A) \7 a ipf3 = ApV Identity II (8) 

Identity II, after using the chain-rule for differentiation, leads directly to the result 

V&{fi=4 IC) VfiS$) ■ (9) 



When this result is combined with Identity I, one easily obtains the equation of motion for 
the vortex density 



dp - V 



(10) 



This continuity equation reflects the fact that the vortex charge is conserved. A key point 
here is that ji is multiplied by the vortex locating 5-function. This means that one can 
replace K in by the part of K which does not vanish as ip — * 0. Thus in the case of a 
nonconserved order parameter one can replace ji K ^ in the continuity equation by 



- 1 r 



TcV 2 ip Ul + rj v 



In the case of a conserved order parameter the current for p is more complicated because of 
the overall gradients acting in K. Because of the standard form of the continuity equation 
Eq . (|10"D , it is clear that one can identify the vortex velocity field as 



/(2) 

Va = ~\ (12) 

where it is assumed that the velocity field is used inside expressions multiplied by the 
vortex locating ^-function. This is the primary result of this paper. It gives one an explicit 
expression for the vortex velocity field in terms of the original order parameter field. In the 
particular case of n = d = 2 one has the more explicit result 



Notice that the result given by Eq. fllCf ) does not depend on the details of the TDGL model, 
only that the equation of motion is first order in time. 

As an important application of the result Eq.flT^D for v consider the velocity probability 
distribution function defined by 

n P(v ) = (n5(v - v)} (14) 

where v is a reference velocity, n = 5(ip)\V\ is the unsigned defect density, and n = 
(n). We determine P using the gaussian closure method |6|-|9| which has been successful 
in determining the scaling function for the order parameter correlation function. The first 
step is to express the order parameter in terms of an auxiliary field in which is assumed, to 
a first approximation, to have a gaussian distribution. In the theory developed in Ref. J7|, 
the relationship between the order parameter and the auxiliary field is given as a solution 
to the classical interface equation 

V 2 J{m) = V'(M)i? (15) 

where the auxilliary field serves as the coordinate labelling the distance to the defect nearest 
to space point f at time t. The solution of this equation for a charge one vortex is of the 
form 

vjj(m) = A(\rh\)m (16) 
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where A(|m|) vanishes linearly with m for small m with the next term of 0(m 3 ). It is then 
easy to show that one can replace ip by fh in the expression for v. One can determine P{vq) 
by first evaluating the more general probability distribution 



G(£, b) = (6(m)6({; - V,m u )5{b - V 2 m)) 



(17) 



since 



where 



with 



and 



n P(v ) = / d n bl[d^\V($\5(v -v(b,$)G(tb) 



v(b,0) 



J {2) (b,Q 

no 



V(£) = — e e pv\pv* p"n 



(18) 



(19) 



(20) 



[n 
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(21) 



In this last expression it has been assumed that the quench is to zero temperature 



so 



that the noise can be set to zero. The gaussian average determining G(£,b) is relatively 
straightforward to evaluate with the result: 



G{t,b) 



(27rS ) n/2 (27rS 4 )™/ 2 (2tt5( 2 ))™ 2 /2 
where, So = ^(m 2 ) is proportional to L 2 , 

5( 2 ) = i-((Vm) 2 ) 



exp 



25( 2 ) 



u\2 



fl,U 



(22) 



(23) 



and 



Sa 



-((V 2 m) 2 ) 
n 



n 



S (2) 



Sn 



(24) 



The quantities So, S^ 2 \ S 4 are determined from the theory for the order parameter correlation 
function. Using this result for G(£,b) in the expression for the probability distribution one 
can use the usual integral representation for the 5-function to perform the integration over 

— * 

the b field to obtain 



n P(v ) = J n d C ( 2 Jg(l)jlv2 ea;p 



2S( 2 ) 



u\2 



(47T 2 S 7) n/2 V&tM 



exp 



(25) 



where 7 = S 4 (rc) 2 , and the matrix M is given by 



M, 



1 



It is straightforward to obtain the rather clean results 

1 



det(M) 



and 



(26) 



(27) 



(28) 



so that 



n P(vo) = J II ^ (27r g( 2))w2/2 exp [-^y E(g) 2 



(47r 2 So7) n / 2 



exp 



"2Z E 



(29) 



The remaining integrals look formidable but can be carried out if one makes a transformation 
from f£ to x„ via 



(30) 



such that 



1 E(0 2 + x v " 



S( 2 ) 



7 ",/3 



(31) 



One easily finds that 
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N aP = Vm 8*+ [ — - l]v%v§ (32) 
L ^/l + (f )7f 2 J 



where 



^2 _ _. /C r(2} _ rr„\2 ^4 



7/^ (2) = (rc) 2 -^ . (33) 

After changing variables from £ to x one obtains 

fS^\ n ' 2 f 1 W 2 J 



V27T,So/ V27TU 2 / [1 + (iTo) 2 /^ 2 ]^ 2 )/ 2 
where J is the remaining dimensionless integral over the x' s which can be be evaluated 
directly as a separable product of gaussian integrals with the simple result J — n\. Since 
P{vq) is normalized to one, we find on integration over v the result 

n °-\^sJ 2«/ 2 r(f + 1) (35) 

which agrees with the result found by Liu and Mazenko || using a more indirect method. 
After using this result for no one finally obtains the result given by Eq. ([[]). This result 
basically says that the probability of finding a large velocity decreases with time. However, 
since this distribution falls off only as v (n+2 ' ) for large vq only the first moment beyond the 
normalization integral exists. This seems to imply the existance of a source of large velocities. 
It seems likely that this is associated with vortex-antivortex final annihilation. 



The determination of So, S^ 2 \ S4 and v requires a theory for the auxiliary field correlation 
function 

C (i2) = -(m(l)-m(2)> . (36) 
n 

There are two theories available and both are of the gaussian closure type assumed above. 
One, due to Ohta, Jasnow, and Kawasaki(OJK) ||12|| , essentially postulates that Cq(12) is a 
gaussian 

C (12) = S e-' 2/{2L ' 2) (37) 
where r = f[ — r 2 . In this case one easily finds that 



v 2 = |W (38) 

where the coefficient of L 2 /t is undetermined in the theory of OJK. Using the theory devel- 
oped in Ref. |1| for n = 2 one finds self-consistently |14| that 



7T 



(Tc) 



2 



e2 = (1 + ^^ < 39) 



where fi = 0.53721.. is the eigenvalue determined within the theory ||13|| . 

One can go forward and extend these ideas to treat two-point velocity correlation func- 
tions and string-like defects (n = d — 1) as will be discussed elsewhere. 
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